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Abstract 

The partial separability of multipartite qubit density matrixes is 
strictly defined. We give a reduction way from N-partite qubit den- 
sity matrixes to bipartite qubit density matrixes, and prove a neces- 
sary condition that a N-partite qubit density matrix to be partially 
separable is its reduced density matrix to satisfy PPT condition. 

PACC numbers: 03.67.Mn; 03.65.Ud; 03.67.Hk 

Recently, an important task in modern quantum mechanics and quantum 
information is to find the criteria of separability of density matrixes. The 
first important result is the well-known positive partial transposition (PPT, 
Peres-Horodecki) criteria[l,2] for 2x2 and 2x3 systems. There are many 
studies about the criteria of separability for the multipartite systems, see 
[3-8]. 

Generally, the common so-called 'separability', in fact, is the full-separability. 
For multipartite systems the problems are more complex, there yet is other 
concept of separability weaker than full-separability, i.e. the 'partial separability', 
e.g. the A-BC-separability, B-AC-separability for a tripartite qubit pure- 
state pabc[&], etc.. Related to Bell-type inequalities and some criteria of 
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partial separability of multipartite systems, etc., see [9-12]. However, we yet 
need to stricter define the concept of partial separability and find the simpler 
criteria. In this paper, first we discuss how to define strictly the concept of 
the partial separability corresponding to a partition. Next, we give a new 
way that an arbitrary N-partite (N> 3) qubit density matrix always can 
be reduced in one step through to a bipartite qubit density matrix . Thus, 
we prove an effective criterion: A necessary condition of a N-partite qubit 
density matrix to be partially separable with respect to a partition is that 
the corresponding reduced bipartite qubit density matrix is separable, i.e. it 
satisfies the PPT condition. Some examples are given. 

Suppose that Pi^-i N is a density matrix for N-partite qubit Hilbert space 
H = ®^ =1 H S , of which the standard basis is {®fLi | i s >} (is = 0, 1). Let 
Z N be the integer set {1, 2, • • • , N} . If two subsets (r) p = {r l5 • • • , r P } and 
(s) N _ P = {si, • • • , s N _ P } in Z N obey 

1 ^ ri < • • • < r P < N, 1 < S! < • • • < s N _ P < N 
(r) P U(s) N _ p = Z N , (r) P n (s) N _ p = 0(1 < P < N) (1) 

where P is an integer, 1^ P ^ N—l, the set {(V)p i ( s )at-p} fo rm s a partition 
of Zjv, in the following we simply call it a 'partition', and for the sake of stress 
we denote it by symbol (r) p || (s) N _ p . A partition (r) p || (s) N _ p corresponds 

to a permutation S {r)pUa)lf _ p = ^ ['[ ^ P + ['[ ^ ^ , by 
which a new matrix P(r) P \\(s) N _ P from p ili2 ... iN is defined now, whose entries 



are 



P(r) P \\(s) { 



jl-'-jN, k\---k N ^jr 1 ---jrpjs 1 ---js N _p, fcrj • • • k rp k sl ■ ■ ■ k s N _ p (^) 



For instance, Pa\\bcd = Pab\\cd = Pabc\\d = Pabcd, and [pAc\\BD\ ijkhrstu = 
[pABCD] ikjlirtsu , [pc\\ABD\ ijkl rstu = [pabcd] kij i,. t rsw etc " Generally, p (r)p \\ Wjv p ^ 
p ili2 ...i N , unless (r) p \\ (s) N _ p just maintains the natural order of Z N (i.e. 
(r) P = (1, • • • , P) , (s) N _ p = (P+ 1, • • • , iV)), then p(r) P \\(s) N _ P = Phi 2 -i N - 

Lemma. For any partition (r) p || (s) N _ p , P(r) P \\(s) N _ P is still a N-partite 
qubit density matrix. 

Proof. We only consider the case of tripartite qubit, the general cases 
are completely similar (also see [11]). Notice the permutation Sb\\ac , then 
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we have 



Pb\\ac = SpabcS\ S 





1 



1 

10 
10 



(3) 



S is an unitary matrix, therefore Pb\\ac is still a tripartite qubit density 
matrix. □ 

Now, we consider how to more strictly define the partial separability. Ob- 
viously, if a partition (r) p || (s) N _ p maintains the natural order of Z N (i.e. 
(r) P = (1, 2, • • • , P) , (s) N _ P = (P + 1, P + 2, • • • , TV)), then p Wp ||(s )jv _ P = 
Phi 2 -i N un der the standard basis {<8>^Li | i 3 >}, now the (r) p — (s) N _ p - 
separability can naturally be defined as that if Pi x i 2 ...i N can be decomposed 
as P(r) P \\(s) N _p = Phh-i N = J2PaPa,(r) P ® Pa,(s) N _ P with probabilities p a , 

where p a ,{r) P and p a ,(s) N _ P , respectively, are a P-partite and a (N — P)- 
partite qubit density matrixes acting upon ®^ =1 P m and ®^Z^H n for all 

to be (r) p — (s) Ar _ p -separable. However, if the nat- 

(s) N _ P ( i.e. Si < r P ), then 



a, then we call p ili2 ... iN 



ural order of Ziv has been broken in (r) 



generally P(r) p \\(s) N _ p 7^ Phi 2 -i N > the case I s different from the above. For in- 
stance, we consider a normalized pure-state pabcd — | ^abcd ^abcd |j 
I ^abcd >£ Pa <8> Pb <E> P<? <S> Pd of four spin-| particles A, B, C and 
D. Now, assume that | ^> abcd > has a special form as | ^> abcd >— 
52 c ik£ji I U > ® | Jb > <8> I k c > <S> | /d >, where Q fe ,c j7 G C 1 . 

i,j,k,l=0,l 

If we keep up to use the original standard basis, then we cannot directly 
see the partial separability, because this choice of basis is unsuitable. If 
we choose other nature basis {| %a > <8> | kc > <8> \ Jb > ® \ Id >} ( this, in 
fact, means that we are using Pac\\bd), under which we can consider the state 
I ^'acbd >=l ^ac > <E> | $bd > , where | * A c >= 52 c ik \ i A > ® \ k c >, 

i,fc=0,l 

I >= 52 c ji I Jb > ® I h > • Now, pac\\bd = Pac ® Pbd, where , 

j',«=o,i 

Pac =| ^ac >< ^ac \,Pbd =| *bd >< |- | ^ abcd > and I ^ABCD >, 
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in fact, are the same in physics, therefore to call pabcd AC-BD-separable 
is completely reasonable. Similarly, for the rest. Generalize to the cases of 
mixed-states, thus we can generally define the concept of partial separability 
as follows. 

Definition. For the partition (r) p || (s) N _ p , a N-partite qubit density 
matrix pi r i 2 ...i N acting upon H= <S>^ =1 H S is called to be (r) p — (s) N _ p -separable 
if the corresponding density matrix P(r) p \\(s) N _ p can be decomposed as 

P(r) P \\{s) N _ P = J2 Papa '( r )p ® P^Mn-p ( 4 ) 

a 

where p a ,(r) P an d Pa,(s) N _ P , respectively, are a P -partite and a (N — P)- 
partite qubit density matrixes acting upon ^ l=1 H rm and ®n=i Hs n for all 
a, andO <p a <l, J2p a = 1- If Phi 2 -i N ^ n °t ( r ) p — ( s ) N _ P - separable, then 

a 

we call it {r) p — (s) N _ p -inseparable. 

For the distinct partitions pi r i 2 ...i N can have distinct separability. Of 
course, if a Pi x i 2 ...i N is partially inseparable for some partition, then it must 
be entangled. Here, in passing, we point out that how to find the gen- 
eral relations between the partial separability and the ordinary separabil- 
ity (full-separability), generally, is not a simple problem. For instance, 

we can make such a multipartite qubit density matrix P (similar to the 
theorem 1 in [13,14]), and by using of the technique in this paper, we 

can prove that p always is partially separable for all possible partitions 
{r) p || (s) N _ P (1 ^ P ^ N — 1) , but P is entangled (not full-separability). 

In order to find the criteria of partial separability, first we discuss how to 
reduce a multipartite qubit density matrix in one step through to a bipartite 
qubit density matrix. For a given partition (r) p \\ (s) N _ p , let two sets (r) p 
and (s) N _ p , respectively be separated again as follows, 

(r')p, = {r[, • ■ ■ , r' P ,} , (r") p „ = {r'{ , • • • , r" p „ } , one of them can be the null set 

(s')o' = {s[, • ■ • , Sq/} , (s")q„ = {s'l, • • • , Sq»} , one of them can be the null set 

r i < r 2 < ' ' ' < r P'-> r l < r 2 < ■ ■ ■ < r P „ (5) 

s[ < s' 2 < ■ ■ ■ < s'q, , s'{ < s' 2 ' < • ■ ■ < s'q,, 

(r) p = (r') p , U (r'V , W) p , fl (r") pl , = 0(0 < P', P" < P and P' + P" = P) 

{s) N _p = (s') QI U (s") Q „ , (s') Q , n (s")q„ = (0 < Q', Q" < N - P and Q' + Q' = N - P) 
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now we rewrite the partition added these partitions as [(r') p , , (r") pil ] 
Now we define the matrix or, n ,,,, 1 ,,r, n ,, M 1 by 

^[{r ') P l,( r )p"JII[(« )m-P'A s )m-P"\ J 



P [{r') P „{r") P ,,]\\[{s>) m _ P „{s")m- P »] 

with form as [p] 3 



Jxix 2 —x N , yiy2—DN 



the submatrix in p^...^ consisting of all entries 

(6) 



which must be a 4x4 matrix, where the values of x k and (k — 1, • • • , N) , 
respectively, are determined by 



x k = i for k G {r') p , , x k = 1 — i for k G {r") pri 

x k = j for k G (s')q' , %k = 1 - j for k G (s") Q „ 

y k = u for k G (r') P , , y fe = 1 - u for fc G (r") p „ 

= v for k G (s')q' , Uk = 1 - ^ for k G (s")q» 



(7) 



where -u,i> = 0, 1. E.g. 



form as [p\ iji(1 



P[{ac),0]\\1(b),(d)] = the submatrix in pabcd consisting of all entries with 

MoooijOooi Hoooi.oioo Hoooi.ion [p1oooi,iiio 
Moioo,oooi Hoioo.oioo Hoioo.ion [p]oioo,iiio 
[p\ 1011,0001 Ip] 1011,0100 Ip] 1011,1011 l/ 5 ] 1011,1110 
[p\ 1110,0001 Ip] 1110,0100 Ip] 1110,1011 l/ 5 ] 



iji(l—j),uvu(l—v) 



(8) 



J 1110,0100 LrJ 1110,1011 Lrj 1110,1110 J 

etc.. Now we define the 4x4 matrix p/,^ ,^ ^ by 



P {{r) P -{s) N _ P ) 



and P\u, 



P {{r) P -{s) N _ p ) 
for all possible [{r>) P , ,{r") P „] \\[{s') Q , ,{s") Ql , 



P [(r>) P ,,{r") P »]\\[{s') m _ P ,,(s") 



[(r')p/,(r")p»]ll[(^) tn _ P /,(.") r , 



-p" 



are not repeated 



(9) 

where we notice that there are indeed repeated P[( r0p/ ,(,») p „]||[( s0m _ p; ,( s » )m _ p// ]; 

m fact, P[ ( r') p „(r") P „}\\[(s') m _ p/ ,(s") m _ pl/ } 
- P[{r>) P ,,(r") P „] II [(«")„_ W«') ra -P'] 

= P[(r") p/ 'Xr') p/ ]\\[(s') rn _ P ,,(s") m - P "] = P [(^")p»,(r')p']ll[(^") m -P»,(^) m -P']' 6tC - F ° r 

instance, we have 



P(A-BC) 



P[(A)M\\[{BC)M + P[(A),0]\\l(B),(C)] 
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P(B-ACD) = P[(B)M\\[(ACD)$\ + P[(B) ,$]\\UC) ,(D)] + P[(B),0]||[(AD),(C)] + P[(B),0]||[(A),(CD)] 
P(AC-BD) = P[(AC),9]\\[(BD)M + P[{AC) M\\[{B) ,{D)\ + P[(A),(C)]\\[(BD),9)] + P[(A) ,(C)]\\[(B) ,{D)] 
P(AC-BDE) = P[(AC)M\l(BDE),9] + P[(AC),9]\\[(BD),(E)] + P[(AC),$]\\[(BE),(D)] 

+Pl(AC),9]\\[(B),(DE)] + P[(A),(C)]\\l(BDE),0] + P[(A),(C)]\\[(BD),(E)] (10) 

+Pl(A),(C)]\\l(BEUD)} + P[(A),(C)]||[(B),pE)] 

etc.. 

As an example, the above reduction procedures from Pabcd to P(ac-bd) 
can be described as Pabcd — ► P(ac-bd) = P[(ac),<H]\\1{bd),$]+P[(ac),0]\\ K b ), (d)] + 

Pl(A),(C)]\\l(BD)M + P[(A),(C)]\\[(B),(D)] 

= a A +a x +a +a A , where the submatrixes cr A , cr x , a and a A , respectively, 
consist of the entries 'A', 'x ','<>' and 'A' in pabcd as in the following figure 
(o" x is just the matrix in Eq.(8)) 



0000 0001 0010 0011 0100 0101 0110 0111 1000 1001 1010 1011 1100 1101 1110 mi 



0000 A 
0001 

0010 
0011 
0100 

0101 A 
0110 

0111 
1000 

1001 

1010 A 

1011 
1100 
1101 
1110 

mi A 



X 



X 



X 



X 



o 



o 
o 



o 



A 



A 



A 



A 



X 



X 



X 



A 



A 



A 



A 



A 



A 



A 



A 



o 



O 
O 



O 



O 



o 
o 



A 



A 



A 



A 



A 



A 



A 



A 



x 



x 



X 



X 



A 



A 



A 



A 



o 



o 
o 



o 



Similarly, we can consider higher dimensional cases. As for the ordinary 
bipartite qubit density matrix pab, we can take p(A-B) = Pab- 

Sum up, generally we can define the 4x4 matrix P^ r ) P -(s) N P ) ^ or a 
given (r) p || (s) N _ p . In addition, it is easily verified that for any partition 



x 



X 



X 



X 



A 



A 



A 



A 
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(r) P || (s) N _ P , P(( s)jv _ p _( r)p ) is the transposition of P( (m)p _ Wjv p ), therefore 
from viewpoint of partial separability, we don't have to distinguish between 
the partitions (r) p || (s) N _ p and (s) N _ p \\ (r) p . 

Theorem 1. For any partition {r) p \\ (s) N _ p , P(( r ) P -(s) N P ) ^ a bipartite 
qubit density matrix, therefore p^ ^ y in fact, is a reduction of the N- 
partite qubit density matrix Pi x i^..% N - 

Proof. The fact must proved only is that P(( r ) P -(s) N P ) ^ s surely a bi- 
partite qubit density matrix. Here we only discuss in detail the cases of 
quadripartite qubit states, since the generalization is completely straightfor- 
ward. In the first place, we prove that the theorem holds for a pure-state 
Pabcd- Suppose that pabcd —\ * ' abcd >< ^ abcd | is a normalized pure- 
state, where | ^abcd >= J2i :j ,k,i=o,i c ijki I U > <S> \ 3b > <S> \ k c > ® \l D >, 



^2i,j,k,i=o,i \°ijki\ — 1- Let 

I $A >= C W I ix > ® I h I $x >= C *3i0--3) \ i x>®\jy> (12) 

i,j=0,l »J=0,1 

I $o >= c ij(l-i)i I *x > ® I jy >, I ^A >= ^2 ^'(I-Oia-J) I > ® I > 

i,j=0,l ij'=0,l 

where x and y are form particles. Make normalization, we obtain pa —\ 

(fA >< <PA I, | V?A >= ^ I $ A >, Px =| <Px >< <p x |, <Px = ^ I $ x >, 
Po =| <p<> >< </?<> |, | <p >= p" 1 | $ >, p A =| (f A >< ip A |, | <p A >= vX 1 | $ A 

>, where the normalization factors are 



va - j icijiji 2 , Vx — y~] \ c iji(i-j)\ 

i,3=0,l V i,j=0,l 



*7o = |cij(i-i)j| 2 ) »7a = / 5^ I c u(i-»)(i-j)| 2 ( 13 ) 

y »,j=o,i y ij=o,i 

It can be directly verified that from Eq.(lO) we have 

P(AC-BD) = VaPa + vlpx+ vlpo + vlpA (14) 

where pa, Px,Po,Pa & U are bipartite qubit pure-states. It is easily seen that 
since | ^ abcd > is normalized, VA + vl+vl + ill= E \ c ijki\ 2 = 1- This 

i,j',fe,/=0,l 
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means that P(ac-bd) is a bipartite qubit density matrix(a mixed state) for 
this pur-state pabcd- 

Secondly, if Pabcd = Y.PaPa(ABCD) is a mixed-state, where every p a {ABCD) 

a 

is a quadripartite qubit pure-state with probabilities p a , then from Eq.(lO) we 

have P(ac~bd) = E a P«, (p<x)(ac-bd) ■ Since ever y (p<*)(ac-bd) is a bipartite 
qubit density matrix, P(ac-bd) is a density matrix (a mixed-state). 

A similar way can be extended to higher dimensional case, the key is that 
when p iu ... tiN is a pure-state, then P[ {r ') P „(r») p// }\\[(s') m _ P „(s>% 



Jm-P" 



■ ^[{r')p / ,(r")p"]\\[W)r n -P'^") m -P"} >< ^[(r') P ,,(rnP<<]\\l(s') m -pn(sn m -P> / y Where 

the pure-state 

I *[M l -,(r") J -,]||[(y) m _p„(yo m -P»] >= E |*i>®---® |o*> 

ij=0,l 

(xi, rr 2 , • • • ,xn are determined by Eq.(7)) (15) 
therefore we just have 

I *u,-^> (16) 

Yl I *[(r')i-,(r")p«]ll[(-') m -P'.(*") m -P»] > 

for all possible [(r') pl ,(r") pll ] || [(s') Q , ,(a") «] , 
and l*r, n , lure /■> /■//■> 1 > arc not re P eat ed 

By using of this relation, make the similar states as in Eq.(12), and make 
generalization to mixes-states, we can prove that generally, a mixed-state 
density matrix Pi x ...i N can be reduced through to the bipartite qubit density 
matrix P( (r)p _ Ww _ p ). □ 

The following theorem is the main result in this paper, it is an application 
of PPT condition for multipartite qubit systems. 

Theorem 2 (Criterion). For a given partition (r) p \\ {s) N _ p , a nec- 
essary condition of a N-partite(N^ 3) qubit density matrix p ili2 ... iN to be 
( r )p — (s) N _ p -separable is that the reduced bipartite qubit density matrix is 
separable, i.e. P(( r ) p -( S ) N p ) satisfies the PPT condition. 

Proof. We only discuss in detail the case of quadripartite qubit, it can be 
straightforwardly generalized to the case of arbitrary N-partite qubit. In the 
first place, we prove that this theorem holds for a quadripartite qubit pure- 
state. Suppose that the pure-state pabcd is AC-BD-separable. This means 
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that if we choose the natural basis {\ tA > ® \ jc > ® \ r B > ® \ sd >} , then 

Pac\\bd = Pac <8> Pbd, where p A c = | * ac >< * ac |, I ^ac >= J2 c ij I 

«J=o,i 

U > ® | JC >, E l C ij| 2 = !> and =1 >< *BD I) I *BD > = 

M=0,1 

XI ^rs I r B > ® I s d >! E Mrs| 2 = 1- From the above ways, it easily 

r,s=0,l r,s=0,l 

checked that the bipartite qubit density matrix P(ac-bd), in fact, can be 
rewritten as 



P(AC-BD) = CT& + &x + &o + &A = &(AC) <S> <T(BD) + &(AC) <S> & 



BD 



^U)®^ ^*^)®*^) (17) 



where a^ A c) =| $(ac*) >< $(ac) I , we already write | <3>(ac) >= Yl e i 

i=0,l 

i x >,6i = Qj and | > <8> \ i c > — >\ i x >■ Similarly, a f =| ® f >< 



I' I $ f ^ >= £ £ I jx >! £ = ^"(w) and I > ® | (l - i) D >- 



AC AC 



j=0,l 

j'a; >, and similarly for <J(bd), a f B v \ > e ^ c -- Now, P(ac-bd) can be written as 



P(AC-BD) = V(AC)V(BD)P(AC) ® P(BD) + ^(AC)^ v\ P(AC) <8> P/ v\ (18) 



where p ( ac) = (?7(ac)) 1 I $(ac) >< $(ac) |,^(ac) = JE N| 2 - Now > P(AC) 



i=0,l 

is a density matrix of a single particle. Similarly, for P^vy P(bd), P^ i 
Since 

V(AC)V(BD) + tfAQV}^ + tf A vW(BD) + 'f/v^/v 



^(AC) + ^vX V(BD) + V, vX ) = 1 (19) 



therefore P(ac-bd) is a separable bipartite qubit mixed-state. The PPT 
condition for separability of 2x2 systems is sufficient and necessary [2], thus 
P(ac-bd) satisfies the PPT condition. Similarly, for other partial separability. 

9 



Secondly, we prove that this theorem holds yet for partially separable 
mixed-states. Suppose that pabcd is a AC-BD-separable mixed-state, then 
under the same natural basis there is a decomposition as Pab\\cd = ^2 a PaPa(AC) 
Pa(BD) -i where p a (AC) an d Pa( bd) both are bipartite qubit pure-states as in the 
above for all a, < p a < 1, J2p a = 1. From the above reduction operation, 

a 

obviously we have 



P(AC-BD) = ^2Pa [pa(AC) ® Pa(BD)] {AC _ BD) 



(20) 



According to the above mention, every \_p a (AC) ® Pa(BD)] (ac-bd) * s a se P ara " 
ble bipartite qubit mixed-state, this leads to that the convex sum P(ac-bd) 
in Eq. (20) still is a separable bipartite qubit mixed-state, and it must satisfy 
the PPT condition. 

Similarly, we cane prove higher dimensional cases. □ 
Corollary. If the reduced bipartite qubit density matrix (piii 2 ---jiv)((r-) 
violates the PPT condition for a partition (r) p || (s) N _ p , then p ili2 ... iN 



( s )n-i 

is 



\ r )p ~ {s) N _ P -inseparable and entangled. 



It, in fact, is the inverse-negative proposition of Theorem 2. 
Examples. Consider two tripartite qubit density matrixes 



Pabc 







l-x 
4 



Pabc 



l-x 
4 



l-x- 
4 



x 
2 





X X 

' 2 2 







l-x 
4 





-f o 



l-x 
4 





l -f o 



l-x 
4 



l-x 
4 



(21) 
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then we have 

(p'abc)(a-bc) = ( Pabc){b-ac) = Pw (22) 

where pw is the Werner state[l,15] which consists of a singlet fraction x and 
a random fraction (1 — x), 

[Pw] ijyrs = xS ijtrs + J (l - x) S ir 5 js 

<Soi,oi — <Sio,io — ~~ <Soi,io — ~~ 'S'io.oi — 2 (23) 
and all the other components of S vanish. 

It is known [1] that when | < x < 1 pw violates the PPT condition, it leads 
to that p' ABC is A-BC-inseparable and p" ABC is B-AC-inseparable. 

By using of the above theorems and corollary, in some special cases we can 
make a iV-partite qubit from 2 N ~ 2 bipartite qubit density matrixes, which 
is partially inseparable for a given partition. As in the above, for the case 
of tripartite qubit we take two bipartite qubit density matrixes a^, <7( 2 ) and 
real numbers p±, p 2 , < pi, p 2 < 1 such that a = Pi(T(i)+ p 2 C(2) is a bipartite 
qubit entangled state ( then it violates the PPT condition). If we want to 
construct a tripartite qubit entangled state Pabc which is B-AC-inseparable, 
then we can take the entries of Pabc by 

[pABc] ijktrst = Pi hl)] iM r ' fOT k = 1 alld 1 = r 

[pABc} ijkyrst = P2 [<7(2)] jijSr , for k = 1 - % and t = 1 - r (24) 
[pABc] ijkirst = 0, for the rest (i, j, k, r,s,t = 0, 1) 

It can be verified that Pabc is a tripartite qubit density matrix, and is B-AC- 
inseparable. In fact, (pabc)(b-ac) = T w hich violates the PPT condition. 
Similarly, for A-BC and C-AB. The above way can be generalized to obtain 
a ( r )p — (s) 7V _ P -inseparable density matrix from a bipartite qubit entangled 

2 N-2 

state in form as r = Yl > where all are some bipartite qubit density 

i=i 

matrixes. 
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